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CONVOLUTIONS INDUCED DISCRETE PROBABILITY 
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Abstract. This paper proposes another constant that can be associ- 
ated with Fibonacci sequence. In this work, we look at the probabil- 
ity distributions generated by the linear convolution of Fibonacci se- 
quence with itself, and the linear convolution of symmetrized Fibonacci 
sequence with itself. We observe that for a distribution generated by 
the linear convolution of the standard Fibonacci sequence with itself, 
the variance converges to 8.4721359. . . . Also, for a distribution gener- 
Cr . ated by the linear convolution of symmetrized Fibonacci sequences, the 

Q, variance converges in an average sense to 17.1942. . . , which is approxi- 

mately twice that we get with common Fibonacci sequence. 

1. Introduction 

Fibonacci sequence is known to have some fascinating properties. The 
most classical of these is the convergence of the ratio J "ii to the golden 
mean (p PQ,[2]. Viswanath [3], proved that with certain randomness intro- 
duced in Fibonacci sequence, the n th root of the absolute value of the n th 
term in the sequence converges, with a probability of 1 (i.e. with extremely 
rare exceptions, almost surely), to another constant 1.13198824 .... In this 
paper, we propose yet another constant that can be associated with the Fi- 
bonacci sequence. This constant, corresponds to the limit variance of the 
distribution generated by the convolution of Fibonacci sequence with itself. 
The variance of such a distribution asymptotically approaches to 8.4721359 

a: 

The paper is organized as follows: 

Section [2] of the paper discusses, in brief, some of the known constants as- 
sociated with the Fibonacci sequence. In section [3] we propose the use of 
Fibonacci sequence convolution to generate discrete probability distribu- 
tions. In section [J] we discuss the results followed by a few observations in 
section with conclusions in section [6J 

2. Some Constants Associated with Fibonacci Sequence 

It is well known that the classical Fibonacci sequence is generated by the 
recursion given in eqQ] below [4 J 

/[n] = /[n-l] + /[n-2] where/[0] = 0;/[l] = 1. (1) 
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A closed form expression for generating Fibonacci sequence is given by the 
Binet's formula, as in Eq|2l below, that relates golden mean 99 to the se- 
quence. 

L J y/E y/5 

where tp = + 2 ■ This is one of the most important and well known con- 
stants that has been associated with the Fibonacci sequence ever since Ke- 
pler showed that the golden ratio is the limit of the ratios of successive 
terms of the Fibonacci sequence. A decade ago, Viswanath [3] computed 
another constant that explained how fast the random Fibonacci sequences 
grow. In his work, he introduced randomness in the recurrence that gener- 
ates Fibonacci sequence. Introducing randomness in Fibonacci sequence, 
we get a random Fibonacci sequence that is defined by the recurrence 
f[n] = f[n— 1] ± f[n — 2] with signs chosen as given below in Eq. [3J 

,r -j = j f[n - 1] + f[n - 2], with probability 0.5; 
1 f[n — 1] — f[n — 2], with probability 0.5. 

Viswanath showed that in the case of the above random Fibonacci sequence, 



Vl/HI -»■ 1-13198824 . . . with probability 1, as n -*■ 00. (4) 

In this work, we propose yet another constant that is closely related to the 
Fibonacci sequence. 

3. Convolution of Fibonacci Sequence and Discrete 
Probability Distributions 

The motivation for this work comes from attempts to exploit the discrete 
nature of integer sequences for generating discrete probability distributions. 
In [5 J , the authors have looked at slow growing sequences and their convolu- 
tions, that approximates a Gaussian distribution with a mean squared error 
of about 10 -8 or even less. In this section, we look at the discrete probability 
distributions generated by a single convolution of Fibonacci sequences. 
Let xi[n] and X2[n] be two discrete sequences. Then the linear convolution 
of the two sequences is defined as, in Eq. ([5]), below. 



y\n 



y] Xi[k]x2[n-k] (5) 



If xi[n] and X2[n] are two finite length sequences of length L and M respec- 
tively, then the length of y[n] is L + M — 1. If the two sequences to be 
convolved, have the same length L, then the length of y[n] is 2L — 1. In this 
work, both x\\p\ and X2[n] are finite length Fibonacci sequences of the same 
length, L. The convolution result is taken as the profile of the discrete prob- 
ability distribution. The set of indices, n, namely 1, 2, 3, . . . ,L + M — 1, of 
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y[n] = x[n] * x[n], is considered the set of values that a discrete random 
variable X can take. The probability of X = n is defined by 



P(X 



n) 



y\n 



sum(y[n]) 



where, sum(y[n]) 



L+M-l 



(6) 



n=l 



Figure Q] illustrates the convolution profiles obtained for different lengths 
of the classical Fibonacci sequence. Figure [2] gives the convolution of two 
Fibonacci sequences of the same length L, with symmetry employed at L/2. 
In this case, we generate Fibonacci sequence of length L/2 and then sym- 
metrically extend it to length L. Figure [2] also compares the convolution 
result with an estimated normal probability density function with the same 
variance as the convolution result. 




Figure 1 . Linear Convolution of Fibonacci Sequences 



4. Results 

In the previous section, we looked at the linear convolution of Fibonacci 
sequence and the discrete distributions generated. We find that, increasing 
the length of Fibonacci Sequence and employing symmetry results in a profile 
that is similar to a very narrow Gaussian profile. For these distributions, we 
plot the variance and standard deviation as functions of the length of the 
sequence. From figure El we observe that the variance of the distribution 
smoothly converge to a constant. This is due to the asymptotic and rapid 
growth of the Fibonacci sequence for large values of n. We also find from 
figure El that the convergence, in case of Fibonacci sequence with symmetry 
employed, is only in the average sense. 
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Figure 2. Fibonacci sequence (symmetry employed) convo- 
lution and narrow Gaussian distribution 
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Figure 3 . Variation of Variance with Length of Fibonacci Sequence 



5. Observations 

We give below, some of the interesting properties, associated with Fi- 
bonacci sequence, that we have observed. Though, we have established a 
few of them mathematically, we are currently working on the analysis of the 
rest of them. 

• Observation 1. With the standard Fibonacci sequence itself be- 
ing considered as a distribution function, the variance converges to 
4.23606797750108. 
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Figure 4. Plot of Index of Absolute Maximum for convolu- 
tion of Fibonacci Sequences 



• Observation 2. When we convolve (linear) standard Fibonacci se- 
quence with itself and take the resulting sequence as a distribution 
function, the variance converges to 8.47213595500216. 

• Observation 3. The value mentioned in observation (2) is twice the 
value mentioned in observation (1). 

• Observation 4. Let S\ be the standard Fibonacci Sequence. Let 52 
be the sequence S\ in the reverse order. Let S3 be the linear convolu- 
tion of Si and S2. The variance of S3 converges to 8.47213595500216. 

• Observation 5. Let Si be the standard Fibonacci Sequence. Let S2 
be the sequence Si in the reverse order. Let S3 be the linear con- 
volution of Si and S2. Let S4 be a sequence obtained by linearly 
convolving S3 with Si or S2. The variance in either cases saturates 
to 12.7081989582623. This value is 3 times the value listed in obser- 
vation (1). 

• Observation 6: The maximum value in the sequence resulting from 
the linear convolution of two standard Fibonacci sequences of length 
L occurs at 2L — 2. 

• Observation 7: Let Si be a standard Fibonacci, L length sequence. 
Let S2 be the sequence that is symmetrically extended version of 
Si. S2 has a length 2L. Convolving S2 with itself yields another 
sequence S3, of length 4L — 1. Considering S3 as a distribution, the 
variance of S3 converges in an average sense to 17.19423665579735. 
The swing is between 17.4442399455347 and 16.9442333660600. 

From figure [U we find that the sequence, resulting from the linear con- 
volution of two increasing Fibonacci sequences of the same length, has its 
maximum at the extreme. From figure [H we find that the maximum is 
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located at 2L — 2. Mathematically, one can prove this by many ways. How- 
ever, we follow the approach of proof by contradiction. For this, it is enough 
to show that 

(I) y[2L - 3] < y[2L - 2], (II) y[2L - 1] < y[2L - 2] 

where y[2L — 1] is the last element of the sequence y[n] of length 2L — 1. This 
is because the linear convolution of two monotonically increasing functions 
will have only one absolute maximum. 

Proof. First we prove y[2L — 1] < y[2L — 2]: 

y[2L - 1] = {f[L])\ (7) 

=►= (/[L-l] + /[L-2]) 2 (8) 

y[2L-l}= (/[L-l]) 2 + (/[L-2]) 2 + 2/[L-l]/[L-2] (9) 

Now, we look at y[2L — 2]. 

y[2L-2] = f[L-l]f[L] + f[L-l]f[L] (10) 

= 2f[L-l]f[L] (11) 

We need to prove that Eq.([9|) is less than Eq. ffTUI) 
Let us assume that y[2L — 1] > y[2L — 2] i.e., 

(f[L - l]) 2 + (f[L - 2])) 2 + 2f[L - l]f[L - 2] > 2f[L - l]f[L] (12) 

>2f[L-l].(f[L-l]+f[L-2]) (13) 

>2/[L-l] 2 + 2/[L-l]/[L-2] (14) 

(f[L - l]) 2 + (f[L - 2]) 2 + 2f[L - l]f[L - 2] > 2(f[L - l]) 2 + 2f[L - l]f[L - 2] (15) 
From (|15p we find that, 

(/[L-2]) 2 >(/[L-l]) 2 (16) 

This implies that f[L — 2] > f[L — 1] which is a contradiction arising due to 
our assumption that y[2L — 1] > y[2L — 2]. Therefore, 

y[2L- 1] <y[2L-2] VL>3 (17) 

We now look at the first inequality y[2L — 3] < y[2L — 2]. Let us assume that 
y[2L-3]>y[2L-2]i.e., 

y[2L - 3] = 2f[L - 2]f[L] + (f[L - l]) 2 (18) 

=> 2f[L - 2}f[L] + (f[L - l]) 2 > 2f[L - l]f[L] (19) 

Substituting for f[L] = f[L - 1] + f[L - 2], we get, 
2(f[L - 2]) 2 + 2f[L - l]f[L - 2] + (/[£ - l]) 2 > 2(/[L - l]) 2 + 2/[i - l]f[L - 2] (20) 

2(/[L-2]) 2 >(/[L-l]) 2 (21) 

^V5>^| (22) 
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This is a contradiction, as the ratio of \ L ~J approaches the golden ratio 

if which is greater than v2- Thus we find that on convolving Fibonacci 
sequence of length L, (L > 3), with itself, results in a sequence y[n], that 
has its maximum at n = 2L — 2. □ 

From the previous section H] we note that the variance of a discrete distri- 
bution generated by the linear convolution of Fibonacci sequence with itself, 
saturates to a constant of value 8.4721359. . . . For a distribution generated 
by the linear convolution of symmetrized Fibonacci sequences, the variance 
saturates in an average sense to 17.1942. . . , which is approximately two 
times 8.4721359. . . . 

6. Conclusion 

This work proposes another constant 8.4721359. . . that can be associated 
with Fibonacci sequence. The constant corresponds to the value to which 
the variance, of the discrete distribution generated by linear convolution 
of Fibonacci sequences of the same length, converges to. Further, it is 
interesting to observe that on performing linear convolution of a Fibonacci 
sequence of length L with itself, the maximum value always occurs at 2L — 2. 
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